Two-Scale Annihilation 
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The kinetics of single-species annihilation, A + A ^ 0, is investigated in which each particle 
has a fixed velocity which may be either itv with equal probability, and a finite diffusivity. In one 
dimension, the interplay between convection and diffusion leads to a decay of the density which is 
proportional to t~^^^. At long times, the reactants organize into domains of right- and left-moving 
particles, with the typical distance between particles in a single domain growing as t^^*, and the 
distance between domains growing as t. The probability that an arbitrary particle reacts with its n^^ 
neighbor is found to decay as n~^^^ for same- velocity pairs and as n"^''* for H — pairs. These kinetic 
and spatial exponents and their interrelations are obtained by scaling arguments. Our predictions 
are in excellent agreement with numerical simulations. 



Single-species diffusion-controlled annihilation, A + 
A 0, exhibits classical mean-field kinetics when the 
spatial dimension d > 2, in which the concentration 
c{t) decays ast~^, and nonclassical dimension-dependent 
kinetics for d < 2 with a slower concentration decay, 
c{t) oc [1-7]. In one dimension, the geometric re- 

striction to nearest-neighbor interactions leads to rela- 
tively large departure from the mean-field kinetics, as 
well as a spatial organization of reactants. In this well- 
studied case, it is known that c{t) asymptotically de- 
cays as {Dt)~^/'^, independent of the initial concentra- 
tion. The complementary situation of single-species an- 
nihilation where the reactants move ballistically has re- 
cently begun to receive attention [8-12]. Perhaps the sim- 
plest example is the deterministic ± annihilation process, 
where each particle moves at a constant velocity which 
may be either +v or —v [8,9]. When the densities of the 
+v and —V particles are equal, c{t) decays as (co/vt)^/^. 

In this letter, we consider single species annihilation 
when the particle transport is a superposition of convec- 
tion and diffusion — we term this system the stochastic 
± annihilation process (Fig. 1). Although the concen- 
tration decays as when only one of the transport 
mechanisms — either convection or diffusion — is oper- 
ative, the combined transport process leads to a faster 
concentration decay of t^^^^. Our goal is to understand 
this unusual decay law and its attendant consequences 
on the spatial distribution of reactants. While there has 
been fragmentary mention of some aspects of this sys- 
tem [7,10], here we give primarily new results and a self- 
contained account of the basic phenomena. 

To set the stage for our approaches and results in the 
stochastic ± annihilation process, it is first helpful to pro- 
vide a simple derivation for the decay of c{t) in the deter- 
ministic ± process. Let us consider a system where par- 
ticles are placed with concentration Cq in a box of size L, 
and denote by c{L,t) the time dependent concentration. 
Initially, there are N = cqL particles, and the difference 
between the number of right and left moving particles is 
of the order of AA^ = |iV+ - iV_| ~ Vn. Eventually, aU 
particles who belong to the minority-velocity species are 



annihilated and thus c(L,t — oo) ^ AN/L ^ (cq/LY^^. 
We assume a scaling from for the concentration, c(L, t) ~ 
(co/i)^/^/(z) with z ~ L/vt. According to the above ar- 
gument, f{z) — + const, in the z — > limit. Conversely, in 
the short time limit, 2 — *■ c», the concentration cannot 
depend on the box size, so that f{z) must be proportional 
to z^/^. Thus we find 



(1) 



As a consequence, the system organizes into right- and 
left-moving domains whose size is of the order of vt. 

In the diffusive case, either one particle or no particles 
survive the annihilation process in a finite box, depending 
on the parity of the initial number of particles. Follow- 
ing the above line of reasoning, we may write the scaling 
ansatz c{L, t) ~ L~^f{z) with z = L/^/lji. Here the rel- 
evant time dependent length scale is V Dt. In the limit 
z ^ 0, the concentration is independent of L, thereby 
implying f[z) ~ z. Therefore the time dependent con- 
centration is given by 
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The crucial new feature in the stochastic ± annihila- 
tion process is that particles with the same velocity can 
mutually annihilate because of their interaction which is 
driven by diffusion (Fig. 1). A useful way to determine 
the decay in this process is to consider separately the 
role of convection and diffusion on the kinetics. Because 
of the convection, particles organize into right-moving 
and left-moving domains as outlined above. Inside each 
domain, however, diffusive annihilation between same- 
velocity particles takes place. We assume that the diffu- 
sive annihilation mechanism leads to an effective time de- 
pendent "initial" concentration, co(i) ~ {Dt)~^/^, which 
plays the role of cq in Eq. (1). Thus we obtain 
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(3) 
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Intriguingly, the concentration in the stochastic ± annihi- 
lation process is predicted to decay as t~^^^ even though 
c{t) decays as if either diffusion only or convection 

only is the transport mechanism. 
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Fig. 1 Space-time evolution of particles in the stochastic ± 
annihilation process. 

An alternative method to determine the decay law, 
which provides additional insight into the relative effects 
of diffusion and convection, is dimensional analysis. If the 
particle diffusion coefficient is D, then the stochastic ± 
process is fully characterized by the initial concentration 
Co, the velocity v, and D. From these parameters, the 
only variable combinations with the dimensions of con- 
centration are cq, l/vt, and l/\/Di. On physical grounds, 
we anticipate that these three concentration scales should 
enter multiplicatively so that the time-dependent concen- 
tration can be expressed in a conventional scaling form. 
Accordingly, we write the time dependent concentration 
in the form 

in which the dimension of the right-hand side is mani- 
festly a concentration. The exponents p and a can be 
now determined by requiring that the above expression 
for c{t) matches with: (a) the diffusion-controlled behav- 
ior c{t) (Dt)^^^^ for t < ~ -D/w^, which is the char- 
acteristic time below which the drift can be ignored for 
a particle which undergoes biased diffusion; and (b) the 
ballistically-controUed behavior c(t) (co/ut)^/^ when 
t < td — l/(-Dco), which is the time for adjacent parti- 
cles to meet by diffusion. Thus by matching Eq. (4) with 
{Dt)~'^/'^ at r^, one obtains p = 0, and then matching 
Eq. (4) with {cQ/vtfl'^ at td gives a = 1/2. This then 
reproduces Eq. (3). 

To test this decay law, we performed Monte Carlo sim- 
ulations using the following realization of the reaction 
process. Initially all sites are occupied with either a -|- or 
a — particle with equal probabilities. A simulation step 



consists of picking a particle at random and moving it a 
single lattice site in the direction of its velocity. If the 
target site is occupied, then both particles are removed 
from the system. Time is updated by the inverse of the 
number of particles. The simulation was carried up to 
10^ time steps on a periodic chain of 10^ sites and an av- 
erage over 10^ realizations was performed. The data for 
c{t) is strikingly linear over a substantial time range on 
a double logarithmic scale (Fig. 2). The local two-point 
slopes of the data in the time range 10^ ^ t ^ 5 x 10'* 
give an exponent value of 0.745. We interpret the con- 
stancy of this data as evidence that the actual value of 
the exponent is 3/4. It is worth noting that a Fade anal- 
ysis of the exact short-time power series gives an estimate 
for the decay exponent of approximately 0.72 [13]. This 
provides a rough estimate for the magnitude of the varia- 
tion of the effective exponent between the early time and 
asymptotic regimes. 
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Fig. 2 Simulation data for the concentration (•) versus time 
on a double logarithmic scale. A line of slope —3/4 is shown 
for reference. 
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Fig. 3 Definition of the basic distance scales that character- 
ize the spatial organization in the stochastic ± annihilation 
process. 

Having established the decay exponent numerically, it 
is of interest to consider the consequences of this unusual 
decay law on the spatial distribution of reactants. In par- 
ticular, since c{t) decays as t"'^/^, one might expect that 
the average separation between nearest-neighbor parti- 
cles grows as t^/'^. However, if there remains any vestige 
of the domain organization that is associated with the 
deterministic ± process, then more than one length scale 
may be needed to characterize this spatial distribution. 
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Such multiscale behavior has been observed previously in 
diffusive two-species annihilation [14] and the associated 
consequences lead to new insights about the system. To 
investigate possible multiscale behavior in the stochastic 
± annihilation process, we introduce the following dis- 
tance scales (Fig. 3): 



{x-+{t)) ^f'-+, {xaon.it)) r^f^^ 



(5) 



which are defined to be, respectively, the average distance 
between neighboring same-velocity pairs, H — pairs, — h 
pairs, and the average length of a domain of same velocity 
particles. 
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Fig. 4 Simulation data for the basic interparticle distances 
{x++{t)} ~ t^++ (o), {x+-{t)) ~ f-^- (A), {x-+it)) ~ r-+ 
(□), and (a;dom(i)} ~ fd"™ (y). Lines of slopes 3/4 and 1 are 
also shown for reference. 

Our Monte Carlo data for these length scales ex- 
hibits considerable curvature on a double logarithmic 
scale (Fig. 4). Thus to estimate the asymptotic behav- 
ior, we studied the systematic variation of the slopes of 
linear least-squares fits as the data at the earliest times 
are progressively eliminated. The effective exponents ob- 
tained in this manner vary considerably; for example, for 
{x++{t)), the eflFective exponent systematically increases, 
but at a progressively slower rate, from 0.699 to 0.734. 
Together with relatively strong numerical evidence that 
the concentration decays as t~^^^, we conclude that the 
actual value of is 3/4. This accords with the expec- 
tation that {x++{t)) should scale as l/c(t). Similar finite 
time corrections occur in the exponent estimates for the 
remaining length scales defined above. For these cases, 
the effective exponent values are all increasing as short- 
time data arc systematically deleted and it appears that 

, V (_, and fdom are all very close to 1, asymptot- 
ically. That is, the corresponding lengths are governed 
by the ballistic particle motion, but again with consider- 
able finite time corrections. The case of {x+-{t)) is espe- 
cially problematic, as the effective exponent changes from 



approximately 0.80 to 0.93 over the time range covered 
by our simulation. Evidently, more extensive simulation 
would be needed to determine the asymptotic exponent 
values unambiguously by simulation alone. 

A new usehil way to characterize the spatial range of 
bimolccular reactions is the collision probability, P(n), 
defined as the probability that the reaction partner of 
a given particle is its n*^ neighbor. Eventually, every 
particle reacts with some collision partner in one dimen- 
sion and the distribution of the distances between part- 
ners provides a measure of the reaction "efficiency" . In 
the deterministic ± process, for example, this probability 
can be obtained analytically [8,9,15]. Let us denote the 
velocity of the n^^ neighbor by = ±1, and the local 
velocity sum by Sn — X]"=o'^«- ^ right moving particle 
initially at the origin reacts with its (2n + 1) "'-neighbor 
if: (a) S'i > for / = 0,1,..., 2n, and (b) S'2„+i = 0. 
This quantity is precisely the same as the first-passage 
probability for a random walk which starts at the ori- 
gin to return to the origin for the first time after 2n 
steps. Because of this equivalence to an exactly solu- 
ble first-passage problem [16], one has P{2n) = and 
P(2n+1) 2-2"-i(2n)!/n!(n + 1)!. In the limit n ^ oo, 
the probability that a given particle collides with its n*''- 
neighbor is given by 



P(n) oc 



(6) 



Motivated by this power law dependence, we assume, 
in general, that P{n) ~ n~^. The exponent 7 can be 
related to other fundamental exponents of reaction pro- 
cesses, namely, the concentration decay exponent a, de- 
fined by c{t) ~ t~°', and the correlation exponent (3, de- 
fined by £_{t) ^ fP. Here ^{t) refers to the distance over 
"information" about the reactants spread. In a time t, 
only particles within a domain of linear size ^(t) are eli- 
gible to react and thus the surviving fraction, or concen- 
tration, is 



c{t) 



/ dnP{n) 
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(7) 



Consequently, we find the exponent relation 

7=1 + a//3. (8) 

For the deterministic ± process, a = 1/2 and (3=1 
[8,9], and the exact 7 = 3/2 of Eq. (6) is recovered. 
As an illustration, consider, for example, single-species 
diffusion-limited annihilation. The decay and correlation 
exponents are a = 1/2 and (3 = 1/2, leading to 7 = 2 
from Eq. (8). Preliminary simulations appear to confirm 
this result. Similarly, for two-species annihilation, a is 
now equal to 1/4 while (3 remains 1/2 so that 7 = 3/2. 

Let us now consider the behavior of the collision prob- 
ability in the stochastic ± annihilation process. In this 
case, the existence of two length scales in the system sug- 
gests that it is necessary to make a distinction between 
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reaction events that involve particles of the same and of 
different velocities. We therefore define a ballistic cor- 
relation scale ^+-{t) ^ t^*^ , with /?+_ = 1, which is 
associated with H — collisions, i.e., annihilation events 
between opposite velocity particles. Invoking the scal- 
ing relation Eq. (8), we thus find {I) ~ with 

7-1 = 7/4. Similarly, there is a diffusive length scale 

~ t^++ , with = f/2, corresponding to anni- 
hilation events between same-velocity particles. In this 
case, Eq. (8) gives 7++ = 5/2. To summarize, we obtain 

P{n) ~ P+-{n) ~ P++(n) = P—{n) ~ n'^'^. 

(9) 

This behavior is consistent with our Monte Carlo simu- 
lation data (Fig. 5). Notice that over large distances, an- 
nihilation between opposite velocity particles dominates, 
as one would naively expect. 




Fig. 5 Simulation data for P-| — (n) (o) and P++(n) (□) on a 
double logarithmic scale. Lines of slopes —7/4 and —5/2 are 
drawn as guides to the eye. 

Our results can also be generalized to arbitrary spatial 

dimension > 1. In this case, it is necessary to ascribe a 
finite, non-zero radius R to the particles so that there is a 
finite collision cross section for particles to actually meet. 
Let us consider the anisotropic system in which particles 
undergo isotropic Brownian motion with diffusivity _D, 
and a drift along the x axis only, with the velocity taking 
on the value -izvx with equal probability. In t he ballistic 
limit {D = 0), the concentration decays as ^ c^)/ R'^~'^vt 
(since the process is quasi-one-dimensional, the t^^/^ de- 
cay of the true one-dimensional system is still obeyed). 
In contrast, for diffusion-controlled annihilation (?; = 0), 
the concentration decays as {Dt)~'^/'^ for d < 2, and as 
{R'^^^Dt)~^ for d > 2 (with logarithmic corrections at 
the critical dimension d = 2) [7]. Repeating the analysis 
detailed previously for the one-dimensional case in the 
derivation of Eq. (3) from Eqs. (1) and (2), we find the 
concentration decay 



r (i?2-2<i „2)-l/4 i-(d+2)/4 ^ < 2; 

c{t)r^ I {RDv)-^/H-' [ln{Dt/R^)]^^\ d = 2; (10) 
[iR^'i-^Dv)-^/H-\ d>2. 

The combined diffusion and ballistic transport does not 
change the mean-field nature of the annihilation kinet- 
ics when d > 2 and the classical decay is recov- 
ered. For sufficiently low spatial dimension, however, the 
non-classical behavior arises in which the decay exponent 
a = {d + 2)/4:. Thus in low spatial dimensions, the in- 
terplay between convection and diffusion provides more 
effective mixing than diffusion or drift alone, and leads to 
a larger decay exponent than ajiff = d/2 and abaii = 1/2 
which arise when only one transport mechanism is oper- 
ative. 

In summary, the stochastic ± single-species annihila- 
tion process exhibits a t^^^^ decay of the concentration. 
This is faster than the t"^/^ decay that arises when only 
one of the constituent transport processes in the stochas- 
tic ± process, either diffusion and deterministic ± con- 
vection, is present. A microscopic understanding of this 
decay law is lacking, and it seems that a technique beyond 
those typically used to solve one-dimensional reactive sys- 
tems would be needed for the stochastic ± annihilation 
process. At long times, the system exhibits a spatial or- 
ganization in which diffusion controls the short distance 
behavior and convection controls the large distance be- 
havior. We have also introduced the concept of the colli- 
sion probability, P{n), the probability that a given parti- 
cle is annihilated by its n*''-neighbor. For the stochastic 
± process, this probability is further discriminated by an- 
nihilation by same-velocity and opposite velocity pairs. 
These two probabilities decay as P++(n) ~ n~^/^ and 



P+_(n) 



respectively. It will be interesting to 



study the collision probability in other reaction processes 
such as diffusive driven single-species annihilation. 
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